A non-extensive approach to the time evolution of Lyapunov coefficients 
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We study sporadic randomness by means of a non-extensive form of Lyapunov coefficient. We 
recover from a different perspective the same conclusion as that of an earlier work, namely, that the 
ordinary Pesin theorem applies (P.Gaspard and X.-J. Wang, Proc. Natl. Acad. Sci. USA 85, 4591 
(1988) ). However, our theoretical analysis allows us to organize the numerical calculations so as to 
reveal the slow transition from a temporary form of non-extensive thermodynamics, corresponding 
to the prediction of a recent paper ( M. Buiatti, P. Grigolini, A. Montagnini, Phys. Rev. Lett 
82, 3383 (1999)), to the ordinary extensive thermodynamics. We show that the transition takes 
place with a slow decay corresponding to the regression from a non-equilibrium initial condition to 
equilibrium condition. 
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I. INTRODUCTION 

The complex and erratic behavior of chaotic trajecto- 
ries is closely related to the so called sensitive depen- 
dence on initial conditions. The degree of this sensitivity 
is measured by the Lyapunov exponent [1] and this is the 
reason why the Lyapunov exponent is so popular in the 
field of deterministic chaos. The Lyapunov coefficient es- 
tablishes, so to speak, the rate of time increase of the 
function defined by: 



general condition. In fact, according to [3], [4] and [5], 
the conventional form of Eq.(2) has to be replaced by: 



£(i) = [l + (l-Q)A Q t]^, 



(3) 



Ait 

£(t) = lim ; { . 



(1) 



where Xq denotes the generalized (non-extensive) Lya- 
punov coefficient. 

The index Q is referred to as entropic index [2], since 
the authors of [3], [4] and [5] prove the structure of Eq. 
(3) to be a natural consequence of the non-extensive form 
of entropy proposed in 1988 by Tsallis [6], which is ex- 
pressed in fact as a function of Q as follows: 



The meaning of this function is as follows. Wc con- 
sider a given trajectory moving from the initial condition 
x(0), and we refer to it as trajectory of interest. Then we 
study another trajectory, departing from an initial con- 
dition very close but distinct from that of the trajectory 
of interest. We refer to this new trajectory as auxiliary 
trajectory. The quantity Ax(t) denotes the distance be- 
tween the auxiliary trajectory and the trajectory of inter- 
est at time t. Consequently, Aai(0) denotes the distance 
between the two distinct initial conditions. The ordinary 
Lyapunov perspective rests on the assumption that: 



Q-l 



(4) 



£(t) = esp(At), 



(2) 



where A denotes the conventional Lyapunov coefficient. 

The adoption of the non-extensive thermodynamics 
advocated by Tsallis [2], however, shows that the expo- 
nential sensitivity behind the definition of the conven- 
tional Lyapunov coefficient is a special case of a more 



We note that when the entropic index Q gets the value 
Q = 1, the new structure of Eq. (3) becomes equivalent 
to the ordinary structure of Eq. (2) thereby showing that 
the ordinary Lyapunov exponent must be identified with 
Ai. 

The purpose of this paper is to adapt the Lyapunov 
approach to this new perspective with a special atten- 
tion to the case of an intermittent map of the same kind 
as that behind the dynamic approach to Levy processes 
of ref. [7]. The authors of [7] proved that the dynamics 
within the laminar region of the intermittent map ad- 
here to the prescription of Eq. (3) with Q > 1. This 
means that the function £(f) of Eq. (3) diverges at a fi- 
nite time, thereby implying an exit from the laminar into 
the chaotic region at much earlier times. Note that the 
trajectory leaving the laminar condition enters a region 
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characterized by Bcrnouilli randomness, and so by ordi- 
nary Lyapunov coefficients. This means, in other words, 
that the resulting dynamics are a balance between a reg- 
ular motion described by Eq. (3) and a random motion 
corresponding to the prescription of Eq. (2). Although 
this random motion is sporadic, it has the important ef- 
fect of erasing the memory of the initial condition [8]. 
We want to discuss the "thermodynamic" consequence 
of this sporadic randomness: This is the process of mem- 
ory erasure necessary for the resulting diffusion process 
to become equivalent to an ordinary Levy process [8] . 

The outline of this paper is as follows. Section II is 
devoted to illustrating, with the help of heuristic argu- 
ments, two distinct prescriptions to evaluate the Lya- 
punov coefficient in the non-extensive case. In Section 
III we supplement the heuristic choice of the proper form 
of Lyapunov coefficent by using the generalized version 
of the Pesin theorem [9] established in an earlier work 
[10]. In Section IV we shall show that the assumption of 
the existence of a smooth invariant distribution implies 
the long-time dynamics of the intermittent map here un- 
der study, namely, the Maneville map [11], to be exten- 
sive. In Section V we study the transition from an off- 
equilibrium condition to the natural invariant distribu- 
tion: The corresponding time evolution of the Lyapunov 
coefficients is monitored by using a numerical approach. 
The physical significance of the results of this paper is 
illustrated in Section VI. 



II. TWO DISTINCT PRESCRIPTIONS TO 
EVALUATE LYAPUNOV COEFFICIENTS 

Let us now discuss two different approaches to the time 
evolution of Lyapunov coefficient, which turn out to be 
equivalent the one to the other only in the extensive case. 
For simplicity we shall limit our discussion to the one- 
dimensional case, and more specifically, to the case where 
the phase space is given by the interval [0,1]. Let us 
assume that 



N-l 



Xn+1 = $(X„) 



(5) 



We shall see that the two distinct definitions rest on the 
use of only one auxiliary trajectory and of a conveniently 
large number of them, respectively. The first definition, 
therefore, seems to be closer to the spirit of the function 
£(t) Of Eq. (1). 

We discuss first the case of only one auxiliary trajec- 
tory. At any time step we can define the corresponding 
conventional Lyapunov coefficient as 



A„ = In 



Ax(n + 1) 
Ax(n) 



(6) 



Thus we can define also the time dependent Lyapunov 
coefficient: 



A(N,x ) = A «- 



(7) 



71=0 



We note that this Lyapunov coefficient can also be writ- 
ten under the form: 



N-l 



A(N,x ) =lnj[\ $'(av. 



(8) 



n=0 



where $ (x) denotes the derivative of <f>(x) with respect 
to x. We leave explicit the dependence of A (AT, x) on the 
initial condition x on purpose. In fact, in this paper 
we are interested in the physical effects caused by this 
form of memory, stemming from the laminar region of 
the phase space. As we shall see in Section IV, the phase 
space of the Manncville map is divided into two parts, 
a laminar region, responsible for these memory effects, 
and a chaotic region. A generic trajectory explores the 
chaotic region only from time to time, thereby making 
sporadic the action of randomness. This sporadic action 
produces the interesting effect of erasing the memory of 
the laminar region, and we plan to discuss the "thermo- 
dynamic" consequences of this process of memory era- 
sure. In the simple case where the whole space phase is 
random and mixing, we expect, in agreement with the 
Pesin theorem [12], that the system reaches quickly the 
steady condition: 



IKS 



limjq- 



A(A,x ) _ JqIu | $'(x) | p(x)da 



N 



Jo P( x ) dx 



(9) 



However, there are cases where it takes extremely long 
time for this steady condition to settle. One of these 
cases is given by the Manneville map [11]. This map, 
which is the dynamical system studied in this paper, is 
the prototype of sporadic dynamical processes, and its 
algorithmic complexity has been discussed by Gaspard 
and Wang [13]. 

Let us now illustrate how to define the Lyapunov coef- 
ficient by means of a conveniently large number of auxil- 
iary trajectories. This is the procedure adopted in the pa- 
per by Benettin and Galgani [14]. We run first of all the 
trajectories of interest moving from the initial condition 
x . Then, we have recourse to a first auxiliary trajectory, 
departing from a different initial condition, at a distance 
Ax(0). We move in time by one step. The distance be- 
tween the two trajectories is now Ax(l) > Ax(0). We ar- 
tificially reduce this new distance to the earlier distance, 
denoted by the symbol d. Then we start a new auxiliary 
trajectory, and so on. It is evident that this prescription 
yields as a natural consequence the definition: 



JV-1 



A(A,x ) =Y, h 



Ax(n + 1) 
d 



(10) 



For practical purposes this approach is better than the 
former, since in this case it is not necessary to set Ax(0) 
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arbitrarily small to ensure that at any time step the mod- 
ulus of Ax(n) is, in turn, so small as to safely adopt the 
linearization assumption behind the definition itself of 
the Lyapunov coefficient. Thus, in this case we are nat- 
urally led to 



JV-l 



A(N,x ) = ^ In \ (x n ) 



(11) 



n=0 



As earlier mentioned, the proposal of Eq.(8) rests on 
the use of only one auxiliary trajectory, and the proposal 
of Eq.(ll) involves the use of infinitely many auxiliary 
trajectories. In the ordinary extensive case here under 
discussion these two distinct proposals yield the same 
result thanks to the property: ln(a ■ b) = Ina + Inb. As 
we shall see in a moment, this equivalence is broken in 
the non-extensive case. 

In the case of fractal dynamics [3] the proper entropic 
index Q gets values different from the prescriptions of 
ordinary statistical mechanics. This means that Q / 1. 
We thus find natural to adopt a mobile entropic index 
q, which in principle, can run from q = — oo to q = 
+oo. The necessity of considering values q ^= 1 breaks 
the equivalence between Eq. (11) and Eq. (8). Note 
that the entropy of Eq. (4) implies the conventional Inx 
to be replaced by ln q x = - 1 _ . Thus, we see that the 
^-generalization of Eq. (8) is: 



A q n \N,x ) 



1 



1-q 



'N-l 



i 



s n=0 



• (12) 



We see also that the ^-generalization of Eq. (11) is: 
rjv-i 



Kf\N,x ) 



1 



q 



E (J -i)) 



71=0 



• (13) 



It is evident that the two definitions are equivalent only 
in the case q = l. 

At first sight one would be tempted to believe that the 
definition of Eq.(13) is more convenient that the other. 
In fact, assuming ergodicity, in the limit N — > oo this 

definition would make the ratio A\ (N, xq)/N become 
identical to the average over the invariant distribution of 
the following form of local Lyapunov coefficient 



A (s) 
q 



J | $>) p(x)dx 



I p(x)dx 
o 



(14) 



However, if this definition were adopted, in the time 
asymptotic limit the time dependence of A, '(N, Xo) 
would be linear regardless of the value of q. For this 
reason it is not convenient to use the proposal of Eq. 
(13). In Section III we shall see that Eq. (12) is the 
correct proposal. 



Let us get a preliminary acquaintance with the defini- 
tion of Eq. (12). Let us note, first of all, that 



A( n )(iV,x ) 



-(q-l)A(N,x ) 



(q-i) 



(15) 



This can be easily obtained by evaluating the logarithm 
of the product appearing in Eq. (12): In fact this allows 
us to establish immediately a connection with the time 
dependent Lyapunov coefficient defined by Eqs. (8) and 
(11). For statistical purposes it is convenient to properly 
average the expression of Eq. (15) over the invariant 
distribution. On an intuitive basis one would expect that 
the proper entropy time evolution is given by 



H* q (N) ^<AW(N,x) 
with the averaging defined by: 



< ...> 



q = J dxp{x) q 



(16) 



(17) 



The problems to settle are two. First of all we have to 
give a stronger support to the heuristic definition of gen- 
eralized Lyapunov coefficient of Eq. (15). Then we have 
to establish a criterion, based on dynamics, to decide 
which value to assign to the mobile index q. All this will 
be explained in Section III. 



III. GENERALIZATION OF THE PESIN 
THEOREM 

To properly appreciate the importance of the expres- 
sion of Eq. (16) (supplemented by Eq. (17)) the 
reader should first of all recall the importance of the 
Kolmogorov-Sinai (KS) entropy [15,16]. The authors of 
Refs. [3], [4] and [5] have already pointed out the im- 
portance of expressing the KS entropy in terms of the 
Tsallis entropy [6] , rather than in terms of the Gibbs en- 
tropy. However, they did not go through a derivation 
entirely based on a single trajectory. The generalization 
of the KS entropy within the single-trajectory spirit of 
the Kolmogorov approach involves the evaluation of the 
following form of entropy: 

H q (N) = 1 -^---^- 1 ^"^) g > (18) 

where p(u ...ujn-i) is the probability of finding the cylin- 
der corresponding to the sequence of symbols wo...Un-i 
[18]. This sequence is generated by dynamics in the fol- 
lowing sense. First of all, the phase space is divided into 
cells, each of them is assigned a label lj t ., then we run 
a virtually infinitely extended trajectory, with M time 
steps, and with M » N. This means that the trajec- 
tory is much longer that the time range explored by the 
subsequent entropic analysis. At any time step this tra- 
jectory occupies a given cell, with a given label, thereby 
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generating a virtually infinite sequence of symbols. This 
"infinite" sequence is explored with a moving window of 
size N, and for any window position a string wo---Wjv-i is 
detected. We have to establish the frequency of occurence 
of this string, to determine thus the corresponding prob- 
ability p(wq...ujn-\), which is then used to fix the value 
of the entropy of Eq. (18). This is the reason why we 
say that this entropy is generated by dynamics. 

The KS entropy is an entropy per unit of time. This 
means that the proper generalization of the KS entropy 
should be given by 



h q = UrriN- 



MP. 

' N ■ 



(19) 



We are thus in the right position to address the problem 
of definining in a non-ambiguous way the earlier men- 
tioned dynamic entropic index Q. If we chose an ar- 
bitrary q the existence of the limit of Eq. (19) would 
not be guaranteed. In the case where we set q = 1 the 
generalized definition of KS entropy becomes identical to 
the ordinary definition, and the KS entropy is related to 
the positive Lyapunov coefficients of the trajectory by 
means of the Pesin theorem [9]. If this limit does not 
exist, namely, it is either a vanishing or an infinite quan- 
tity, one would be tempted to conclude that the dynamic 
process under study is incompatible with thermodynam- 
ics in the sense of Kolmogorov. However, if it happens 
that a given value of q exists, such that the limit of Eq. 
(19) exists and is finite, then we consider it to be the 
dynamic entropic index of the process under study. We 
call Q this "magic" value of q, which is in fact the proper 
dynamic entropic index of the process. In an earlier work 
[17] a numerical calculation was made to establish Q in 
the case of a text of only two symbols, with strong cor- 
relation. The case of many more symbols is beyond the 
range of the current generation of computers. 

This arduous task is made less difficult by the fact that 
recently a proper generalization of the Pesin theorem has 
been established [10]. According to the prescription of 
this paper the entropy H q (t) of Eq.(18) reads 



H g (t) 



1 - 51- 1 J dxp{x)mt 1 x) 1 -i 



1 



(20) 



where the symbol t denotes time regarded as a continuous 
variable. In fact, when the condition TV >> 1 applies, it 
is legitimate to identify N with t. The function p(x) 
denotes the equilibrium distribution density and S the 
size of the partition cells: According to [10] the phase 
space, a one-dimensional interval, has been divided into 
W = 1/5 cells of equal size. 

We have to remark that the reason why the Kol- 
mogorov entropy is so popular is because it affords a way 
of establishing the randomness of a process without any 
arbitrary dependence on the way the observation process 
is carried out. In fact, in the case where the partition 
of the phase space into cell is generating [18], the result- 
ing entropy rate is independent of the partition adopted. 



This important aspect is retained by the generalized ex- 
pression of Eq.(20). In fact, this equation apparently 
suggests that the independence of the entropy rate of the 
partition adopted is guaranteed only when q = Q = 1. 
Actually, a recent work [19] has assessed that if the in- 
variant distribution is multifractal the expression of Eq. 
(20) becomes independent of the size of the cells, and 
the proper dynamic entropic index can be different from 
unity. On the other hand, this means that in the case 
where the invariant distribution is smooth, rather than 
multifractal, the dynamic entropic index Q must become 
equal to the unity. We shall discuss the physical conse- 
quences of this important conclusion in the next sections. 



IV. THE MANNEVILLE MAP 

The Manneville map [11] reads: 

x n+ i = $(x n ) =x n + x z n (modl) (z > 1). (21) 

It is characterized by a laminar region given by the in- 
terval [0, d] with d defined by the equation: 



d z +d = l. 



(22) 



It was recently shown [7] that in this case the prediction 
of Eq. (3) holds true with 



Q = l + (z-l)/z. 



(23) 



Note that the condition z = 1 , making the Manneville 
map equivalent to the Bernoulli shift map [1], yields, as 
it must, Q = 1 and thus a form of ordinary statistical 
mechanics. It is easy, as well as convenient, to write 
explicitly the two Lyapunov coefficients of Eqs. (12) and 
(13) in this case. Due to the fact that = 2 wc 

obtain that the two prescriptions result in 



Af)(A) 



and 



9(1-9) 1 
A^(N) = N- — - 



(24) 



(25) 



respectively. We see that Eq. (24) results in a linear time 
evolution of the Lyapunov coefficent only when q = Q = 
1. For q > Q the Lyapunov coefficent increase is slower 
and for q < Q is faster. A steady entropy increase per 
unit of time can only be defined for q = Q = 1. We make 
therefore the conjecture that this property migth hold 
true even when Q ^ 1, in agreement with the result of 
numerical work of [17]. We also see, as already pointed 
out in Section II, that Eq. (13) yields Eq. (25) and 
this time dependent Lyapunov coefficent turns out to be 
linear in N regardless of the value of q used. 

To help the reader to understand the significance of 
the results of this paper, it is convenient to review the 
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physical arguments used in [7] to make the prediction of 
Eq. (23). The authors of [7] used an intermittent map 
with two laminar regions, each of which is equivalent to 
the laminar region of the Manneville map. The sojourn in 
one of the two laminar regions means uniform motion in 
the laboratory frame of reference with vclociy W in cither 
the positive or the negative direction. The random walker 
created by this dynamic process makes jumps of intensity 
| a; | in either the positive or the negative direction with a 
probability II(|a;|), which is related to the waiting time 
distribution ip(t) in one of the two equivalent laminar 
regions by means of the key relation 



n(M) = 



w 



(26) 



On the other hand the intermittent map used in [7] yields 
a waiting time distribution ip(t) whose explicit expression 
is: 



V>(i) =d z - l [l + d z - l {z-l)t]- 



(27) 



We will focus our attention on the interval 1.5 < z < 2 
of the map parameter. This is the range of z for which 
only the first moment of ip(t) is finite and anomalous 
diffusion is expected [7]. For 1 < z < 1.5 also the sec- 
ond moment of ip(t) is finite and ordinary diffusion is 
obtained. Finally, in the case of z > 2 all the moments 
are infinite, thereby conflicting with the requirement of 
producing diffusion compatible with an Hamiltonian de- 
scription [20]. 

The maximization of the non-extensive entropy 
^(ndxl)), namely, the non-extensive entropy expressed 
as a functional of FI(|a;|) yields an analytical expression 
for II(|a;|) which is compatible with the inverse power law 
form inherited by n(|a;|) through Eq. (26). This is the 
origin of the prediction of Eq. (23). The authors of Ref. 
[7], using dynamic arguments, prove that the Manneville 
map yields 



C(t) = [i-(z-iK- 1 t]- 



(28) 



As far as the power index is concerned, this dynamical 
property is compatible with Eq. (3) and with the predic- 
tion of Eq. (23). Unfortunately, the coefficient Aq of Eq. 
(3) through comparison with Eq. (28) turns out to be 



Aq = zx 



z-l 



(29) 



This dependence on the initial condition conflicts with 
the fact that the waiting time distribution is independent 
of the initial condition, being a statistical average based 
on the assumption of a random but uniform injection 
from the chaotic to the laminar region [21]. 

The rationale for this conflict is that the generalized 
Pesin theorem of Section III implies that the dynamic 
entropic index gets the ordinary value Q = 1. In fact, it 
is well known [13] that the invariant distribution is given 

by 



p(x) = 



const 



x- 



-l 



(30) 



This means a smooth distribution, and consequently ac- 
cording to the arguments of Section III, Q = 1. This 
prediction seems to be in a harsh conflict with the pre- 
diction of Eq. (23). However, as we shall see in the next 
section, it is not so, since the prediction of Eq. (23) refers 
to a microscopic process with high memory, before the 
crossing from the laminar to the chaotic region occurrs. 



V. NUMERICAL RESULTS 

The prediction Q = 1 of the generalized version of 
the Pesin theorem of Eq. (20) refers to a condition of 
equilibrium with the system in its invariant distribution. 
To establish a contact with the prediction of Eq. (23) it 
is convenient to average the Lyapunov coefficient of Eq. 
(12) on a non-equilibrium distribution density. Thus we 
use the numerical work to evaluate 



A q (N) =<A q n \N,x) > r 



(31) 



Note that this prescription has not to be confused with 
that of Eq. (16) even if both are prescriptions to make an 
average over the same Lyapunov coefficient, that of Eq. 
(15). Here, < ... > ne denotes an average different from 
that of Eq. (17), since we use p{x) rather than p q {x). 
Furthermore the statistical weight p(x) is not the natu- 
ral invariant distribution of Eq. (17). We select on the 
contrary, as a statistical weight, a uniform distribution 
from i = 0toj; = A<d. All the numerical calcula- 
tions of this section refer to A = 10~ 4 . Note that, as 
discussed in [7], the range of ^-values corresponding to 
the emergence of Levy processes is given by the interval 
[1.5,2]. We select the intermediate value z = 1.7 and, 
for computational convenience, the Lyapunov coefficent 
of Eq. (15) is expressed in the equivalent form 



l-^(N,x)^-^ 



1 



(32) 



We aim at establishing at any given time region the value 
of the entropic index q which makes the corresponding 
non-extensive entropy increase linearly with N. This is a 
temporary "magic" Q, since , as we shall see, the windows 
of linear increase have a finite time duration. We study 
the time derivative of A q (N) of Eq. (31), namely: 



B q (N) =A q (N + l)-A q (N). 



(33) 



The purpose of the numerical calculation is that of de- 
termining the time regions where B q (N) is constant. 

The theoretical remarks of Section IV imply that only 
one proper "thermodynamical" region exists, this being 
given by B q (N) constant for an infinite time interval. 
This is the region of ordinary statistical mechanics with 
Q = 1. However, we shall show that the time evolution 
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of the Manneville map ensuing a given off-equilibrium 
initial condition goes through non-extensive phases of in- 
creasingly time duration, and this time duration becomes 
larger and larger as the deviations from the extensive con- 
dition become slighter and slighter. 

In Fig. 1 we plot the quantity B\(N). Note that 
B\ (TV) is the difference between two subsequent iteration 
steps of the Lyapunov coefficient A (TV, x) of Eq. (8) av- 
eraged over a non equilibrium distribution. We see that 
a relatively extended time region shows up, ranging from 
TV = to TV ~ 900, where B\(N) is approximately con- 
stant with an almost vanishing value. This is the region 
where according to the rule established intuitively with 
the help of Eq. (24) the existence of a vanishing non- 
extensive Lyapunov coefficient at a given q implies that 
a finite non-extensive Lyapunov coefficient might exists 
at a larger value of q. In fact, as we shall see with the help 
of Fig. 3, at the magic value of q given by Eq. (23), this 
time region is found to correspond to a finite Lyapunov 
coefficient. 

After this initial time region Bi(N) undergoes an 
abrupt increase followed by a slow regression to a con- 
stant non vanishing value. This fits very well the the- 
oretical prediction of Section IV. The slow regression 
corresponds to the time evolution of the initial non- 
equilibrium distribution towards the final invariant dis- 
tribution. This has the effect of realizing a condition 
equivalent to that of Eq. (9). In conclusion, Fig. 1 sheds 
light into the process of transition to equilibrium , and 
gives a further evidence to the fact that this equilibrium 
is of extensive nature. This figure suggests that the final 
equilibrium condition is reached at about TV ~ 15,000. 
However, as we shall see with the help of Figs. 2, 4, 5 and 
6, at this time the process of regression to equilibrium is 
still active. 



_ 0.3 - 




1 5000 20000 



10000 
N 

FIG. 1. Bi(N) as a function of the iteration step TV. Bi(N) 
has a vanishing value until TV ~ 900 after which undergoes 
an abrupt increase followed by a slow regression, with sta- 
tistical fluctuations, to the constant value predicted by Pesin 
theorem. 



We think that a condition of genuine regression to 
equilibrium implies that all the trajectories, departing 



from the selected off-equilibrium initial distribution, have 
crossed the border between the laminar and chaotic re- 
gion at least once. The process of regression to equilib- 
rium is slow because the motion of trajectories close to 
x = is slow. To shed light on this important aspect, let 
us adopt the continuous time approximation to Eq. (21) 
and let us write: 



dx/dt = x z . 



(34) 



In line with the prescription of denoting by t the discrete 
time TV when TV >> 1 applies, from here throughout the 
remainder of this Section, we shall be denoting time with 
symbol t. The solution of Eq.(34) for a trajectory with 
initial condition x(0) is given by 



x(t) = [x(O) 1 - 2 - (z- l)t]Th 



(35) 



Using this solution it is easy to find the time at which 
the first trajectory, that belonging to the right border 
of the initial distribution, exits from the laminar region. 
This time, denoted by T, is given by 



d x ~ z A 

r -r^TKf) , -- 1 i- 



(36) 



Before this time no trajectory can exit from the laminar 
region. To evaluate the population decrease after this 
time, we notice that: 



rfM = M(0)J, 



(37) 



where M(0) is the number of trajectories within the lam- 
inar region at the initial time. Thus, using Eq. (35) we 
obtain 



dM 
~dT 



M(0) 



A [l + d'-^z-l)^' 



(38) 



which makes it possible for us to establish the time evolu- 
tion of the population at time t for t > T. In conclusion, 
we get 



M(t) = M(0) (t < T) 



and 



M(t) = M (0) 



A [l + d*- l {z-l)t] — 



(39) 



(t > T). (40) 



It is interesting to remark that the time derivative of 
TVf(t), at t > T, as resulting from Eq. (40), turns out 
to be proportional to the waiting time distribution ip(t) 
of Eq. (27). This means that M(t), although depend- 
ing on an arbitrary initial condition, at times larger than 
T reflects the stationary and statistical nature of ip(t). 
This fits the numerical observation made herein with the 
help of Fig. 6 that the process of regression to equilib- 
rium is not affected by the return of the trajectories from 
the chaotic to the laminar region. We shall see, in fact, 
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that the long-time behavior of M(t) of Eq. (40) is a fair 
indicator of the process of regression to equilibrium. 

In Fig. 2 we plot M(i)/M(0) as resulting from the 
numerical treatment. The result obtained coincides with 
the theoretical prediction of Eqs. (39) and (40). It looks 
like the mirror image of the curve of Fig. 1, thereby con- 
firming that the relaxation of the non-equilibrium Lya- 
punov coefficient is due to the trajectories crossing at 
least once the border between laminar and chaotic re- 
gion. 




20000 



FIG. 2. Time evolution of the population of the laminar 
region. The numerical result coincides with the time evolution 
predicted by Eqs. (39) and (40) and the initial plateau, lasting 
for the time T ~ 900 coincides with the prediction of Eq. (36). 




15000 



20000 



10000 

t 

FIG. 3. -B1.41, multipied by the factor of 10 3 , as a function 
of time. The value of Q predicted by Eq. (23) turns out to 
be the proper entropic index in the initial time interval [0, T] 
with T ~ 900. The subsequent quick drop to zero signals the 
regression to q = 1. 



In Fig. 3 we study the time derivative of the Lyapunov 
coefficient of Eq. (12) with q = Q = 1.41, that is, the 
value obtained from Eq. (23) with z = 1.7. We see that 
this rate of entropy increase per time unit is constant in 
the initial interval < t < T . This means that the non- 
extensive entropy corresponding to the magic value of the 
entropic index q given by the theoretical prediction of [7] , 



Eq. (23), realizes, as it must do, the Kolmogorov condi- 
tion of an increase linear in time. From Fig. 3 we also 
see that after the first escape from the laminar region, at 
t = T, the rate per unit time of this non-extensive en- 
tropy drops quickly to much smaller values with a subse- 
quent slower relaxation to a vanishing value in the time 
asymptotic limit. This means that in this time region 
the time increase of this non-extensive entropy becomes 
slower than a linear function in time, thereby suggesting, 
on the basis of the intuitive rule established by Eq. (24), 
that the proper entropic index in this case is smaller than 
that given by the prediction of Eq. (23). In fact, we have 
already argued that at equilibrium Q = 1. 




2000 



FIG. 4. Q as a function of time: The fast drop at t = T. We 
see that at t = T, when the first trajectory escapes from the 
laminar region, Q(t) drops from the value 1.41, established 
by Eq. (23), to the value 1.01 in about 200 iteration steps. 
This fast drop is followed by a much slower relaxation to the 
final value of 1, detailed in Fig. 5. 



In Figs. 4 and 5 we show the change of the entropic 
index from the prediction of Eq.(23), Q = l + (z—l)/z , 
corresponding to a trajectory motion not yet affected by 
the sporadic randomness, to the asymptotic value Q = 1, 
produced by the repeated action of the chaotic region of 
the phase space. In section III we have seen that it is the 
smooth nature of the invariant distribution that makes it 
impossible to adopt the general prescription of Eq. (20) 
in a form different from that of the ordinary Pesin theo- 
rem. On the other hand, this invariant distribution is 
realized through the process of escape from the lami- 
nar region, and so by the action of sporadic randomness, 
whose effect is illustrated here in detail. The calculation 
is done by searching, at any time t > 0, for the value of q 
realizing, temporarily, the condition of linear increase of 
the non-extensive entropy. This, in turn, is realized by 
looking for the value of q making B q (t) of Eq. (33) con- 
stant over windows of finite size. The window sizes are 
reported in both figures as small intervals, whose length 
tends to increase as the value of Q decreases. Fig. 4 
shows the remarkable fact that the transition from the 
region where only order exists to that where sporadic 
randomness begins showing up is signalled by the fast 
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drop of the cntropic index Q. 

Nevertheless after the fast drop at t = T, the equilib- 
rium value Q = 1 is reached asymptotically in time with 
a slow regression process. This is illustrated by Fig. 5, 
which shows more clearly than Fig. 4, that the time du- 
ration of this temporary non-extensive thermodynamics 
becomes larger and larger as Q comes closer and closer to 
the equilibrium value Q = 1. We think that the adoption 
of the non-extensive formalism results in an impressively 
sensitive indicator of the slight departure from the final 
equilibrium distribution. 



1.001 | | 

1.0008 - - 
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1.0004 - — 

1.0002 - 

1 4000 8000 12000 16000 
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FIG. 5. Q as a function of time: Details of the slow regres- 
sion following the fast drop. This figure shows, more clearly 
than Fig. 4, the windows of entropy linear increase. The size 
of these windows becomes larger and larger at greater and 
greater times. The largest window shown here corresponds 
to Q = 1.00001. In the scale of this figure smaller values 
of q becomes indistinguishable from the value q — 1, which 
corresponds to a window of infinite time size. 

The process of regression of Q(t) to the equilibrium 
value Q = 1 is closely related to the process of transition 
from the initial unstable distribution to the final invari- 
ant measure. This latter process, in turn, is determined 
by the trajectories exiting the laminar region, and con- 
sequently is related to the time evolution M(t) of Eqs. 
(39) and (40). However, an essential part of this pro- 
cess of relaxation to equilibrium might be played also by 
the trajectories that from the chaotic region are injected 
back into the laminar region. To establish if this true or 
not, it is convenient to monitor numerically the process 
of transition to equilibrium. This is done as follows. The 
interval [0, 1] is divided into C cells of equal size. In our 
calculations we set C = 100. We consider M trajecto- 
ries with initial conditions uniformly distributed over the 
whole interval [0, 1] and we iterate all of them N times. 
We set M = 10, 000 and N = 100, 000. In accordance 
with the prescription of Section III, we denote time with 
the symbol t. At this stage we evaluate how many tra- 
jectories are found in a given cell with the label i. We 
call Mi this number and we set P[ q = Mi/M. This stage 
does not afford yet a proper numerical determination of 
the equilibrium distribution. This is so because with a 



finite number of trajectories M and cells of a finite size 
1/C, the quantity P^ q will turn out to be a fluctuating 
function of the iteration time t = N. The intensity of 
these fluctuations depends on the selected values for the 
numbers C and M. To bypass this limitation we make 
a time average on r further iterations after the time t, 
thereby defining 

<P! q >^\jz^(t); (41) 

t'= o 

where < Pi{t) > T denotes the probability of the z-th 
cell t further iterations after the time t. The result of 
this calculation is scarcely dependent on the value of tau 
adopted if this is much greater than the time scale of the 
fluctuations of P^ q : In our case for t > 500 the right 
hand side of Eq. (41) remains practically constant for all 
the cells. Thus, we can omit the dependence of < P^ q > T 
on t and use Eq. (41), with a given value of r > 500, 
to define our numerical equilibrium distribution, which, 
for the sake of simplicity is again denoted by the symbol 

Now we can address the important issue of the regres- 
sion to equilibrium. First of all, we adopt the same initial 
condition as that used in the earlier calculations. From 
this initial distribution we select a sample of M trajecto- 
ries. At any time step we count how many trajectories are 
found in a given cell, thereby determining Pi(t), namely, 
the probability that a trajectory is found in the cell with 
label i at time t. We compare this probability with the 
equilibrium probability, evaluated according to the ear- 
lier numerical prescription, thereby defining the variable 
Yi(t) as 

Y i (t) = \P i (t)-P? I \ . (42) 

Then we evaluate the relative dispersion Ri{t) of the 
quantities Yi{t) around the equilibrium value P^ q 

Ri(t) = (43) 

i 

We have now to deal with the issue of the fluctuations 
of Ri (t) caused by the adoption of finite values for C and 
M. We follow the same procedure as that adopted to 
determine the equilibrium distribution. This means that 
we make the time average 

< >r= \ E ( 44 ) 

t'= t 

thereby deriving the time average of the relative disper- 
sion Ri(t): 

<Ri(t)>r= <Y fJ q >T - (45) 

i 

We, then, consider a set of numbers R in the interval 
[0, 1] .For each of these numbers we determine the number 
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of iterations T necessary to make the time average of the 
relative dispersion, Eq. (45), smaller than i?,for all the 
C cells. This makes it possible for us to use the function 
R(T) as a fair indicator of the relaxation to equilibrium. 
In fact, for any time T we can say that the corresponding 
distribution departs from equilibrium by an amount of 
the order of R * 100 per cent. The values of r adopted 
range from r = 500 to r = 3000. Within this wide 
interval the change of R(T) is not significant. 

In Fig. 6 we plot R(T) for different values of r. This 
figure shows that at T = 10, 000 and T = 12, 000 the dis- 
tribution stills departs from equilibrium, for all the val- 
ues of r, by a quantity of the order of 3 ~ 4 per cent. We 
need to wait till to a time greater than 20, 000 to detect 
a smaller departure from equilibrium. Consequently, ac- 
cording to the criterion we have adopted, at T = 16, 000 
the departure from equilibrium is expected to be still of 
the order of 3 ~ 4 per cent. 

It would seem to be plausible to make the conjecture 
that the return of the trajectories from the chaotic to 
the laminar region makes the relaxation to equilibrium 
slower than the decay of M(t) given by Eq. (40). This 
is so because we arrived at this expression for M(t) by 
neglecting the process of return of the trajectories to the 
laminar region after the escape into the chaotic region. 
However, the numerical calculation reported in Fig. 6 
shows that it is not so and that M(t) is a good indica- 
tor of the process of relaxation to equilibrium. Thus, we 
can use M{t) to establish if the small departures of Q(t) 
from Q = 1 are a genuine indication of memory of the 
initial condition. For example the long window of linear- 
ity between t w 6, 000 and t w 16, 000, corresponding to 
Q = 1.00001 , is shown by M(t) to correspond to a de- 
parture from equilibrium of the order of 3 ~ 4 per cent. 
This means that this very small departure from the ordi- 
nary condition Q — 1 signals that about 3^4 per cent 
of the trajectories are still in the laminar region. 



0.6 









x = 500 







x = 1000 







X = 1500 


X 




x = 2000 


+ 




x i 





5000 10000 15000 20000 

T 

FIG. 6. R as a function of the time T. Here the number of 
cells C is 100 and the trajectories used ,M, are 10000. 



VI. CONCLUDING REMARKS 

This paper illustrates the dynamical process of the 
memory erasure necessary to realize the processes of Levy 
diffusion, and therefore makes it possible to complete the 
program of the work of [7]. In fact, the authors of that 
paper focused their attention on the microscopic pro- 
cess responsible for anomalous diffusion of Levy kind. 
This implies elementary jumps of a given length x with 
a probability proportional to 1/lxl* 771 . This requires the 
maximization of the non-extensive entropy [7] with the 
entropic index given by Eq.(23). However, when an at- 
tempt is made at establishing a connection between this 
extensive property and the KS entropy, either in the nor- 
mal or in the non-ordinary form [10], a conflict emerges 
due to the fact that according to the work of [13] the 
ordinary version of the Pesin theorem applies. 

Actually the process of regression to equilibrium is 
shown to be very slow and take place according to the 
theoretical prediction of Eqs. (39) and (40). We do not 
find any numerical evidence supporting the plausible con- 
jecture that the process of regression to equilibrium is 
made slower by the process re- injection of the trajectories 
into the laminar region. In conclusion, as an important 
result of this paper, we show that there is a close connec- 
tion between this slow process of regression and the aging 
of the entropic index: The process of aging corresponds 
to the memory erasure necessary to establish the Levy 
statistics [8]. 

This process of memory erasure, beggining at t = T ~ 
900, is not as fast as it seems to be at a first sight. The 
numerical and theoreticaal arguments of Section V prove 
that it is slow. In other words, a very extended regime 
exists, of balance between the randomness of the chaotic 
region and the order of the laminar region, and it is sig- 
naled by values of the entropic index Q very close to 
Q = l, but not coinciding with this value. The reader 
might find surprising that Q = 1.00001 signals the exis- 
tence of memory, in spite of being so close to the predic- 
tion of ordinary statistical mechanichs, Q = 1. Actually, 
our theoretical and numerical arguments prove that this 
very weak deviation from Q = 1 corresponds to a pro- 
cess of relaxation still active. This conclusion, of course, 
is made possible by the fact that the dynamic process 
under study is understood since relatively long time [13]. 

We see essentially two main benefits emerging from the 
results of this paper. The first is that the non-extensive 
entropy might serve the purpose of detecting a residual 
order in the long-time limit of real time series, where we 
would expect a totally random behavior. The second is 
of conceptual nature and has to do with the meaning it- 
self of the non-extensive thermodynamics. The present 
paper suggests that the non-extensive thermodynamic 
condition is not permanent, and rather refers to non- 
equilibrium process of long time duration. 
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